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^ ! Abstract 

^ ' We consider fluid wetting on a corrugated substrate using effective interfacial 

r^ ' Hamiltonian theory and show that breaking the translational invariance along 

«■ 

O ' the wall can induce an unbending phase transition in addition to unbinding. 

Both first order and second order unbending transitions can occur at and out 

^^ ' of coexistence. Results for systems with short-ranged and long-ranged forces 

(<— ^ ' establish that the unbending critical point is characterised by hyperuniversal 

(N' 

^^ ' scaling behaviour. We show that, at bulk coexistence, the adsorption at 

^^ ' the unbending critical point is a universal multiple of the adsorption for the 

correspondent planar system. 
PACS numbers: 68.45.Gd, 68.35.Rh 
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Recently, the subject of fluid adsorption and wetting on structured (non-planar) and 
heterogeneous substrates has began to receive considerable attention |[l|. This work is not 
only a natural extension of studies of wetting on idealised planar surfaces but it is 
also of more fundamental interest since the broken translational invariance along the wall 
necessarily leads to competition between surface tension and direct molecular effects. Thus, 
we may anticipate that new interesting phenomena (phase transitions, scaling, universality) 
will emerge which do not occur for planar systems. In this letter, we report results of 
extensive numerical calculations, supported by approximate non-perturbative analysis and 
scaling theory, of wetting on a periodic (corrugated) substrate. These reveal that novel first 
and second order transitions can take place, directly related to the inhomogeneity along the 
wall. For long-ranged forces, the phase transition, referred as unbending, only occurs for 
sufficiently large wall corrugations (beyond the range of previously employed perturbative 
methods 0), dependent on the wave vector of the corrugation. In contrast, for short-ranged 
forces, the critical threshold is wave vector independent and rather weak. There are three 
aspects of our work that we emphasize in particular. Firstly, the unbending transition 
precedes a wetting (unbinding) transition occurring at a higher temperature (and at bulk 
two-phase coexistence). For second order unbinding transitions, on which we concentrate, 
the location of the wetting transition is unaffected by wall corrugation. Secondly, the location 
of the unbending line and critical point, as well as the interface structure, only depend on 
the amplitude and period of the wall corrugation function through hyperuniversal scaling 
variables analogous to that encountered in the theory of finite-size effects at bulk critical 
points [0. As a consequence, the unbending critical point is associated with non-trivial 
universal amplitude ratios which relate the adsorptions in the non-planar and correspondent 
planar system. Finally, unbending is directly related to nonlinear bifurcation phenomena 
occurring in dynamical systems, a subject whose mathematical aspects continue to attract 
attention [^. 

To begin, we describe the results of a specific mean-field (MF) model of unbending and 
unbinding which also serves to illustrate important scaling properties which we shall later 



put in a more general context. For simplicity, we assume that the wall has a corrugated 
sinusoidal shape ^(x) =acos(qa;), which breaks the translational invariance in one direction 
only. Following the work of earlier authors Jl],^, we take as our starting point the (reduced) 
standard effective interfacial model 
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restricted to the space of periodic solutions which is sufficient for our description of equilib- 
rium phenomena. Here, S is the surface stiffness, W is the binding potential and l{x) is the 
collective coordinate measuring the height of the interface relative to the mean position of 
the wall whose period L satisfies q = 27r/L. We also restrict ourselves to a MF description 
in which the equilibrium profiles Cu are obtained by minimising Eq. (|I]). The importance of 
fiuctuation effects will be discussed later in the context of scaling theory 0. We start by 
considering systems with short-ranged forces at bulk two-phase coexistence and write 

W{E) = -AT e-^ + p e-^^ ; (2) 

so that both the film thickness i and corrugation amplitude a are measured in units of the 
bulk correlation length. With this potential (and positive f3) the planar system undergoes a 
second-order unbinding transition at AT = T.^ — T = such that the MF interface thickness 
and the transverse correlation length diverge at that critical point as £^ ~ — log(AT) and 
^11 ~ AT ^^, corresponding to standard wetting critical exponents /^s = (log) and z/y = 1 
respectively 0. 

For a 7^ 0, the MF configuration (s) are the solutions of the Euler-Lagrange equation 

sc(x) = iy'(4-^), (3) 

solved subject to periodic boundary conditions and where the prime denotes differentiation 
w.r.t. the argument. This deceptively simple looking nonlinear equation can show multiple 
solutions and bifurcations corresponding to different possible phases for the equilibrium 
interface configuration. Whilst a full analytic solution is not possible, it is straightforward 



to show that the solutions exhibit an important scahng property which allows us to collapse 
results obtained for different periods L=2ti/ q onto a universal surface phase diagram. To 
see this, we introduce the new variables ri = i — ip — i,^ and t = qx so that (^ becomes 

f/ = AT 2 ( e"" - e-2^) + a cos t, (4) 

which is the equation of a forced inverted nonlinear oscillator. Here the overdot corresponds 
to differentiation w.r.t. t whilst the temperature, stiffness and substrate periodicity are 



combined in the rescaled temperature variable AT = AT/qv^2/3S. Consequently, any new 
phase transition induced by the corrugation amplitude a is not affected by the value of the 
wall periodicity q which only acts to rescale the temperature deviation from T^. In Fig. 
1, we show plots of the mean interface thickness io, defined as the average <i{x) >^., as a 
function of AT for various a, obtained by numerically minimising Eq. (|T|). It can be seen 
that, whilst the location of the unbinding transition is unaffected by the wall corrugation, a 
new phase transition occurs for corrugation amplitudes a>ac~2.914 and AT> ATc~2.12. 
The surface phase diagram is shown in Fig. 2 and exhibits the termination of the first- 
order phase boundary at an unbending critical point as well as representative shapes of the 
coexisting interfacial phases at the transition. Again, we emphasize the universal value of 
the critical corrugation amplitudes a^- (which is independent of q) whilst the temperature 
shift from T^ satisfies ATc(q) oc q. 

Before we discuss further scaling properties that emerge from the exact minimization 
of (|I]), we describe an approximate treatment of the model which recovers the unbending 
transition and yields relatively good values for the critical point. To this end, we suppose 
that the interface configuration, and consequently the free-energy, can be parametrized by 
two variables by restricting ourselves to profiles of the form i{x) ~ £o + (1 — ^)4^{^)- Thus, £o 
is the average interface displacement whilst e measures the extent of interfacial corrugation. 
The bounding value e = 1 corresponds to a completely flat configuration whereas e = refers 
to a configuration with identical corrugation to the wall. Substituting this parametrized 
profile shape into the Hamiltonian, Eq. (|ID, and minimising w.r.t. £o, we are led to the 



following approximate expression for the dependence of the free-energy F on the interface 
corrugation parameter e, 

2 F{e) = ^{l-ef - /\f^^^ (5) 



Eq2 ' ' 2 ' ' /o(2ea) 

where Iq denotes the modified Bessel function of zero-order. The two terms on the r.h.s. 
represent the competition between the surface tension and binding potential effects which 
are each minimized separately by e = 1 (fiat interface) and e = (corrugated interface) 
respectively. Plots of F{e) for various a moving along the unbending line are shown in Fig. 
3 and illustrate the possibility of phase coexistence between bent and rather flat states for 
sufficiently large a. The locus of the unbending transition in the surface phase diagram 
obtained in this approximate manner is shown as the dashed line in Fig. 2 and agrees 
reasonably well with the exact numerical result. Note that the solutions will only depend 
on AT and a, as in the exact solution. This method also has a distinct advantage over 
previously adopted perturbative treatments P| (involving an expansion about the planar 
system) which, whilst not without merit, cannot handle the occurrence of distinct branches 
{i.e. a bifurcation) in the free-energy 0. We also note that the location of the unbending 
critical point within this approximate non-perturbative method can be determined with an 
elegant graphical construction 0. 

We consider now the same phenomena for systems with long-ranged (dispersion) forces. 
For this case, we use the binding potential 

AT B 

^w = -1^ + ^ (6) 

which again describes a continuous unbinding transition in the planar system as T -^ T^^, P| . 
For this system, the film thickness and transverse correlation length diverge as l^^ ~ AT ~^ 
and .^11 ~ AT~^'^, corresponding to critical exponents /5s = 1 and v\\ =5/2 respectively [^. 
Turning to the non-planar geometry, we make the judicial change of variables ri = {i — iI')/^tt 
and t = qx which again reduces the Euler-Lagrange equation (|]) to that of a forced inverted 
nonlinear oscillator: 



ri = AT^{ — -—)+acost. (7) 

Once more, the two scaling variables AT = 2 AT/Sq^£^ and a = a/ij^ determine the multi- 
plicity of solutions and hence the surface phase diagram. 

Plots of the mean interface position £o vs. AT for different a obtained from the numerical 
minimization of Eq. (|l|) are, in essence, the same as that shown in Fig. 1 for short-ranged 
forces and, therefore, are not presented here. The numerical values for the scaled variables 
at the unbending critical point are Oc ~ 2.061 and AT ^ 8.66 which imply a wave-vector 
dependence adq) oc q~^/^ and ATc(q) oc q^/^ for the critical corrugation amplitude and 
temperature shift, respectively. 

The MF results described above suggest that the location of the unbending critical point 
can be understood using scaling theory. To this end, we suppose that, in the planar system, 
the excess free-energy per unit area contains a singular contribution Ftt ~ AT ^""^ (with 
as = and —1 for the model potential (||) and @, respectively 0). In the non-planar 
system, we conjecture that the corresponding quantity is described by the scaling function 

ATf "^g = AT 2-"s H^(aAT^s^qAT-^ii) (8) 

where W{x,y) is the scaling function whose variables correspond to the hyperuniversal 
combination of lengthscales a/i.,^ and q^y 0. Since the singularity in the free-energy at the 
unbending critical point occurs for AT 7^0, we are immediately led to the prediction for the 
critical corrugation amplitude and temperature 

ac(q) (X q"^ ; ATc(q) (x q^ (9) 

consistent with our explicit results, provided that for short-ranged forces we interpret /^s/z^y 
as zero and not logarithmic. We believe that the existence of a finite critical threshold even 
in the g — > limit is a surprising finding of our work. These scaling ideas can be extended 
to the interface structure at the unbending critical point where the hyperuniversal nature of 
the scaling variables x and y play an important role. Here, we concentrate on systems with 

6 



long-ranged forces for which /3s 7^ where the definition of universal critical amplitudes is 
more straightforward. We suppose that, in the vicinity of the unbending critical point, the 
mean interface thickness in the non-planar system is described by the scaling law 



^.A(£,qe||) (10) 



where A{x,y) is a universal scaling function. As a consequence, precisely at the unbending 
critical point, the mean film thickness £q is a universal multiple of the corresponding planar 
adsorption (at the same temperature). Thus, we define the universal critical amplitude ratio 

R = j- at a = ac(q), Ar = ATc(q) (11) 

which we have numerically determined as i? ~ 1.321 (independent of q) calculated using 
our MF theory with the binding potential (^). Note that the definition of R is equivalent 
to the ratio of adsorptions in the non-planar and planar systems. Other universal critical 
amplitudes can also be defined. For example, at the unbinding critical point, the shift in 
the mean interface height relative to the planar system satisfies 
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R' ^ ^° ^:^^^^^ (12) 

«c(q) 
with R' also independent of q. The advantage of this definition is that it is also appropriate 
for systems in which /3s = (log). We have numerically determined that R' = 0.640 and 
0.156 for the potentials (0) and (^ respectively. 

To finish our article, we make two pertinent remarks. Firstly, we have established that 
for a > ttc the first order unbending transition also occurs out of the two-phase coexistence 
for sufficiently small bulk ordering field h. The result of our numerical calculations for short- 
ranged forces including an additional hi term in the binding potential are shown in Fig. 4. 
The existence of an unbending line extending out of bulk two-phase coexistence is analogous 
to prewetting at (planar) first-order phase transitions. Secondly, we have established that 
unbending also occurs for first order wetting transitions in non-planar systems although the 
scaling behaviour is less obvious. A section of the surface phase diagram in the (T, h) plane 



thus shows both prewetting and unbending hues. While this first appears similar to prefilling 
on a wedge, there are profound and subtle differences between unbending and prefilling 
relating to the order of these transitions and their relation with wetting. In summary, we have 
shown that for non-planar systems an additional interfacial phase transition is associated 
with unbinding. The critical point of the unbending transition exhibits novel scaling and 
observable universal critical properties. Further work should concentrate on more general 
wall shapes, calculations with more microscopic models and also aim to establish whether 
the values of the universal critical amplitudes presented here are substantially affected by 
including fluctuation effects beyond mean-field level. At present, simulation studies seem 
best equipped to answer this latter question although renormalization group analysis may 
be possible. 
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FIG. 1. Film thickness vs. AT for different values of a from the numerical minimization of Eq. 
(^. From below, a/^2 = 0.00, 0.50, 1.00, 1.50, 1.75, 1.90, 2.0605, 2.10, 2.15, 2.20, 2.25, 2.30. All 
distances are measured in units of the bulk correlation length. 
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FIG. 2. Section of the surface phase diagram at bulk coexistence showing the unbending coex- 
istence hne which finishes at the critical point AT ~ 2.12 and a^ ~ 2.914. The solid line represents 
the results of minimizing Eq. (||). The dashed line is the result of the variational approximate 
solution (see text). The vertical line AT = represents the second order unbinding transition. 
Schematic representation of the interfacial configuration on either side of the unbending line are 
also shown. 
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FIG. 3. Free Energy from the expression (g) for differents values of a at the transition tem- 
peratures (See Fig. 2). From above, a = 5.0, 4.0, 3.5, 3.0, 2.75 and 2.46866. The last value is the 
critical value Cc within the present approximation. 
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FIG. 4. Phase diagram of the unbending transition in a system of short ranged forces for 
different values oia/ yj2 = 2.5, 3.0, 4.0 and 5.0, from numerical minimization of Eq. (|l|) (continuous 
lines, left to right). The circle represents the unbending critical point for ac~ 2.914. The loci of 
critical points obtained from the approximate model is represented as a broken line. The triangle 
locates the critical wetting temperature. 
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